Driven harmonic oscillator as a quantum simulator for open systems 
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We show theoretically how a driven harmonic oscillator can be used as a quantum simulator for the 
non-Markovian damped harmonic oscillator. In the general framework, our results demonstrate the 
possibility to use a closed system as a simulator for open quantum systems. The quantum simulator 
is based on sets of controlled drives of the closed harmonic oscillator with appropriately tailored 
electric field pulses. The non-Markovian dynamics of the damped harmonic oscillator is obtained by 
using the information about the spectral density of the open system when averaging over the drives 
of the closed oscillator. We consider single trapped ions as a specific physical implementation of 
the simulator, and we show how the simulator approach reveals new physical insight into the open 
system dynamics, e.g. the characteristic quantum mechanical non-Markovian oscillatory behavior 
of the energy of the damped oscillator, usually obtained by the non-Lindblad-type master equation, 
can have a simple semiclassical interpretation. 

PACS numbers: 03.67.Lx, 03.65.Yz, 32.80.Pj 



I. INTRODUCTION 

Recent years have witnessed a considerable experimen- 
tal and theoretical progress in quantum information sci- 
ence Q1Q,E|. A major factor for the rapid development 
has been the increasing ability to control and engineer 
the quantum states of single and few particle systems. 
The development of the field is expected to continue in 
fast pace and some applications, e.g. for quantum cryp- 
tography, are already commercially available. One of the 
most active contemporary sub-fields of quantum informa- 
tion science and of nanotechnology is the study of quan- 
tum simulators HSHH00. 

Simulating quantum systems by classical computers re- 
mains a challenging task due to the inherent inability of 
classical systems to incorporate the quantum features in 
an efficient way. The development of quantum simula- 
tors, i.e. a controllable quantum system imitating the 
behavior of other quantum systems of interest, holds the 
promise for a generation of powerful means and devices 
to study quantum systems and provides additional ways 
to gain new insight into their peculiar quantum features. 

A notable aspect of recent proposals for quantum sim- 
ulators is their ability to cross the usual boundaries be- 
tween different fields of physics 0, 0, H|, e -g- a certain 
configuration of cold atoms trapped in optical lattices 
has been recently proposed to simulate superstrings 
These, as well as more "conservative" quantum simula- 
tors 0, 0, are often based on cold trapped atoms and 
gases because of the recently developed abilities to ma- 
nipulate their properties in an extremely precise way. 

In this paper we develop a quantum simulator scheme 
for open quantum systems exploiting the recent progress 
in trapped ion technologies [2j- The motivation to de- 
velop quantum simulators for open systems stems es- 
pecially from their central role in the progress towards 
fault-tolerant quantum information processing. Indeed, 
decoherence, which is a key issue for successful quantum 
gate implementation, is a topic typically studied in terms 



of open quantum systems [To| . 

We focus on a quantum simulator scheme for a damped 
harmonic oscillator which is a paradigmatic model of 
open quantum systems [lOj and has a wide variety of 
applications ranging, e.g., from quantum optics an d 
nuclear physics |l2j ] to chemistry 0] . We show how 
a driven harmonic oscillator can be used to mimic the 
behavior of the damped harmonic oscillator. In other 
words, we demonstrate theoretically the possibility to use 
a closed system as a quantum simulator for open quan- 
tum system. 

Our quantum simulator is based on the controlled driv- 
ing of a closed harmonic oscillator by sets of appropri- 
ately tailored electric fields. We show how the sets of 
drives in different frequency regimes give rise to various 
types of dynamical features of the driven oscillator, and 
how these shed light on the effects that the various parts 
of the environment spectrum have on the open system 
dynamics to be simulated. The total non-Markovian dy- 
namics of the damped harmonic oscillator is obtained by 
using the information about the spectral density of the 
open system when averaging over the sets of drives of the 
closed oscillator. 

A well known example of experimental realization of 
quantum harmonic oscillator is a single trapped ion || 
and we discuss the implementation of the simulator with 
this system. The effects of random electromagnetic fields 
on the tra pped ion dynamics has been recently stud- 
ied 0, In the usual approaches, the system 
dynamics is calculated averaging over known stochastic 
properties of the environment fields. Our approach is 
taken from the opposite direction. We ask if it is possible 
to drive the closed oscillator in a controlled way so that 
the effect mimics the open system dynamics. Our results 
highlight the possibility to simulate artificial, engineered 
reservoirs, by controlled sets of drives. In general, reser- 
voir engineering is emerging as an active contemporary 
research topic for the control and the fundamental study 
of open quantum systems 0, 0, 0, |2(j • 
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The quantum simulator approach, already at a theo- 
retical level, makes it possible to gain new insight into the 
complex non-Markovian dynamics of the damped oscilla- 
tor. The simulator approach allows to identify the origin 
of the different dynamical features of the non-Markovian 
dynamics of the open system and sheds light on the role 
that the phase of the random environmental fields plays 
in the system dynamics. Furthermore we show how the 
implementation of the quantum simulator with a trapped 
ion allows to interpret a typical quantum mechanical be- 
havior, usually described with a non-Lindblad-type mas- 
ter equation, by semiclassical means. We also discuss 
the implications of the quantum simulator approach to- 
wards the continuos measurement scheme interpretation 
of the quantum trajectories for non-Markovian open sys- 
tems [21| . 

The paper is organized as follows. Section [H] intro- 
duces the model of a closed harmonic oscillator driven 
by controlled fields and sets the framework for the quan- 
tum simulator approach presented in Sec. IIIII Here, the 
damped oscillator model to be simulated, is presented, 
and we furthermore discuss the new insight that the sim- 
ulator approach brings to light. Section Hvl shows how the 
quantum simulator can be implemented with a trapped 
ion and the discussion concludes the paper in Sec. [V] 



II. DRIVEN OSCILLATOR DYNAMICS 
A. The basic scheme: pulsed driving 

We begin by considering the dynamics of a closed quan- 
tum harmonic oscillator which is driven by a single fre- 
quency field for a set of runs which we label by the fre- 
quency lo of the field. In each run of the set lo, the 
duration of the applied field changes. In other words, 
a single set of runs with frequency lo consists of drives 
with increasing duration of the pulses mapping the os- 
cillator dynamics as a function of time. The frequency 
of the driving field changes for different sets of runs and 
depends on the type of the environment and spectral den- 
sity to be simulated. 

It is important to emphasize that the dynamics of the 
closed oscillator is not monitored (measured) in a con- 
tinuos way. Rather, the scheme for a single run of the 
set lo can be summarized as follows: (i) Switch on the 
driving field at time t\ = 0. The state of the oscillator 
may change suddenly due to a sudden switch on of the 
driving field, (ii) Driving of the oscillator with external 
field of frequency lo for < t < t 2 . (iii) Switch off of the 
driving field at t = t%. The state of the oscillator may 
change suddenly due to a sudden switch off of the driving 
field. With the set of pulses, with increasing t 2 for the 
members of the set lo, one maps the oscillator dynamics. 
The key ingredient for the simulator, for the system we 
study, is related to the effects of switching on and off of 
the field. 

The Hamiltonian of a quantum harmonic oscillator 



driven by a time dependent periodic force is |22( 

H = hw U f a + + hF(t)(a + (J). (1) 

Here, loq is the frequency of the oscillator, (a) is the 
creation (annihilation) operator of the energy quanta of 
the oscillator, and the periodic driving force is 



F(t) 



Acos(ujt + ip) 
\/2mfiLOo 



(2) 



where A describes the amplitude and lo the oscillation 
frequency of the periodic force, m is the mass of the os- 
cillator, and ip is the phase of the driving field. 



B. Off-resonant driving 

We assume, for simplicity, that the oscillator is initially 
in the ground state. In this case it can be shown that the 
heating function (n) = (a^a), for lo ^ loq, can be written 



(n>(t,w) - |«| s 



[l + COS 2 (LOt) 



-2 cos(cot) cos(woi) — 2 — sin(wt) sin(woi) 

LOq 



LO 



sin 2 (u>t) 



(3) 



where the phase of the field is chosen for convenience 
tp = and a = Aj \J2mhio^. 

Figures n ( a ) j (b), and (c) show examples of the time 
evolutions of (n) for lo/lo — 0.0, 0.1, 0.2. In the fol- 
lowing two subsections we discuss in detail the heating 
function dynamics for constant field drive (lo = 0) and 
in the adiabatic regime (w < wo). In these frequency 
ranges, the dynamical features of the driven oscillator 
heating function allow to gain interesting insight into the 
non-Markovian dynamics of the open system (to be dis- 
cussed in Sec. Illlll . 



1. Constant field 

When driving the system with a constant field lo = 0, 
the heating function displays sinusoidal oscillations, see 
Fig-E( a )- This can be explained in the following way: (i) 
Switching on the driving field at time t\ = 0. A suddenly 
switched on constant force displaces the ground state of 
the oscillator creating a coherent state D(a)\0) = \a) 
(ii) Time evolution between t\ < t < t 2 . The coherent 
state oscillates in the trap according to the free evolution 
(no change in (n)). (iii) Switching off the field at t — 
t 2 . The initial displacement of (i) is reversed [D(— a)]. 
Moreover, the effect of the second displacement depends 
on the phase of the oscillation of the coherent state. 
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The final state of the system at time ti (after the sec- 
ond displacement) is 



aexp(— iu>ot)), 



(4) 



and the heating function becomes 

(n)(t) = | -a + aexp(-iu Q t)\ 2 = 2\a\ 2 (l - cos u t). (5) 

The second displacement brings the coherent state back 
to the ground state after full oscillation periods whereas 
half the period of oscillation of the coherent state during 
ti < t < t 2 causes a total displacement of 2a. The heat- 
ing function consequently oscillates between and 4|a| 2 . 
Setting oj = in the general result of Eq. © matches the 
simple result [Eq. and confirms this interpretation. 

2. Driving with small frequency in adiabatic regime 

Let us consider now the case < u> <C u>q. Since 
the driving frequency is much smaller than the oscillator 
frequency, the oscillator follows adiabatically the changes 
of the force. Figures ^ (b) and (c) display examples of 
drivings in this regime for co/ujq = 0.1, 0.2. 

Also here, the switches on and off of the drive field 
are thought to happen instantaneously compared to the 
oscillator dynamics. The simple scheme of the previous 
subsection is replaced with: (i) Switching on the field at 
time t\ — 0. As before, the switched on field displaces 
the ground state of the oscillator creating a coherent state 
D(a)\0) — \a) (ii) Time evolution between t\ < t < ti. 
The coherent state oscillates in the trap according to the 
free evolution (no change in (n)). Moreover, the periodic 
driving force with frequency lj changes slowly. The os- 
cillator follows adiabatically the consequent slow motion 
of the oscillator center since w c C wo. (iii) Switching off 
the field at t = t-},. Because of the slow oscillation of the 
driving force, the magnitude of the second displacement 
depends now also on driving frequency u> and on time t: 
D{— cos(wi) exp(— ibj t)a). 

Following this scheme, the state of the oscillator after 
the second displacement can be written 



I — iu>ot 



cos(o;£)a), 



(6) 



and the heating function becomes 

(ri)(t,u) = \a\ 2 [l + cos 2 (uit) - 2 cos(ujt) cos{uj t)] . (7) 

This matches the general result of Eq. J3J) when the terms 
containing lo/luq <C 1 are neglected (adiabatic limit) and 
confirms the interpretation presented here. 

C. On resonant driving 

The equation © for the heating function is not well 
defined for lo — luq. Beginning from the Hamiltonian 





FIG. 1: Time evolution of the heating function of the closed 
driven oscillator in three different frequency regimes for four 
different driving frequencies: (a) uj/uiq — 0, (b) oj/ujq = 0.1, 
(c) lu/u)q — 0.2, (d) lu/ojo — 1.0. Please note that we have 
multiplied the result in (d) by a factor of 0.01 (a weight factor 
given by the spectral density we use in Sec. Mill . 



of Eq. Q one can show that on-resonant driving of the 
oscillator gives for the heating function (for convenience 
the phase of the field is set ip — 0) 



(n)(t,uj Q ) 



Ut 2 



■ uiot sm(2uJot) 



sin 



(wot)] 



Figure ^(d) displays an example of the heating function 
dynamics in this case. We note that the oscillatory terms 
play a minor role contrary to the off-resonant constant 
field and adiabatic cases. 



III. DRIVEN OSCILLATOR AS A QUANTUM 
SIMULATOR 

In the previous section we have shown the main dy- 
namical features of the pulsed driven oscillator in various 
frequency ranges of the driving force (constant force, adi- 
abatic regime, on-resonant driving). We show now that 
by combining the dominant dynamical features from the 
three main frequency regimes, one can construct a quan- 
tum simulator for an open quantum system. We begin by 
recalling the basic features of the system to be simulated. 



A. System to be simulated: damped harmonic 
oscillator 

The damped harmonic oscillator is one of the paradig- 
matic models used to describe the dynamics of open 
quantum systems . We briefly recall the form of the 
master equation for the reduced system density matrix 
which allows to describe the non-Markovian dynamics of 
the damped oscillator. For a more detailed presentation, 
see e.g. Ref. [23j and references therein. 
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The dynamics of a harmonic oscillator linearly coupled 
to a quantized reservoir, modelled as an infinite chain of 
quantum harmonic oscillators, is described, in the secu- 
lar approximation, b y m eans of the following generalized 
master equation [23ll24l] 

dp(t) A(t)+ 7 (t) 



dt 2 
A(i)- 7 (t) 



\2ap(t)a) — a) ap{t) — pit)a* a\ 
\2a) p{t)a — aa) pit) — p(t)aa t ] , 



(9) 



with p(t) the reduced density matrix of the system 
harmonic oscillator. The time dependent coefficients 
A(t) and j(t) appearing in the master equation are 
known as diffusion and dissipation coefficients, respec- 
tively mm. 

The diffusion coefficient for high temperature reservoir 
k B T '/u>o S> 1 and to second order in the dimensionless 
coupling constant g, can be written [23l l2flj 



A(t) = 2g 2 k B T 



{1 



1 + r 2 
(1/r) sm(w *)] }■ 



[cos(wo^) 



(10) 



Above, r = uj c /luq is the ratio between the environment 
cut-off frequency lu c and the oscillator frequency u)q, k B 
is the Boltzmann constant, and T is the temperature. 
For the high temperature case the dissipation coefficient 
7 <C A and therefore is neglected. 

A typical environment of open systems is described by 
an Ohmic reservoir spectral density with Lorentz-Drudc 
cut-off H 



J H = — .,2 



(11) 



for a schematic view, see Fig. [5] The spectral distribution 
is given by 



1{uS) = J(ui)[n e (ui) + 1/2] 

,2 

— coth (u/KT), 



^7. 



7T (jjz -f- UJ 



(12) 



where n e (uS) is the population of the environment mode 
of frequency ui, and Eq. lfTT|l has been used. For high T, 
Eq. (fT5|> becomes 



2k B T 



TT UJ 2 + LU 2 



(13) 



The central parameter r = uj c /^>o describes how on- 
resonant the oscillator is with the reservoir. When r > 1, 
intensive part of the environment spectrum overlaps with 
the oscillator frequency and the decay coefficient A(t) > 
for all times. Consequently, the master equation is of 
Lindblad-type. When r < 1 (Fig. the most intense 
part of the environment spectrum lies in small frequency 
range and the on-resonant intensity is small. In this case 
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FIG. 2: Schematic view of the spectral density and the split- 
ting of the spectrum (r = 0.1). Part I is the low frequency 
part which produces the sudden displacements of the oscilla- 
tor state when switching on and off. Part II is the on-resonant 
part which produces Markovian type heating. 



the decay coefficient A(t) acquire temporarily negative 
values and the master equation is of non-Lindblad-type 

m 

The full solution of the the master equation 10 can 
be found, e.g., in Refs. m,m|. For an initial Fock state 
| to = 0) it turns out that, in the secular approximation, 
the state of the system at each time t is a thermal state 
with "effective" temperature changing in time according 
to the behavior of (n(t)) (see Appendix A): 



Pk,k 



<»(*)> 



lit)) 



1 



1 



i(t)) + r 



(14) 



It is worth emphasizing that the short time non- 
Markovian effects arc described by the non-monotonic 
changes in (n) [see the exact result for (n) in Fig. 3 (d)]. 
The non-Markovian regime is then followed by Marko- 
vian linear heating and finally for long times the system 
reaches a steady state in thermal equilibrium with its 
environment (see Ref. [28] for a discussion of the origin 
of the oscillations of the heating function in the non- 
Markovian regime and Ref. 29] for a discussion of the 
Markovian and thermalization regimes). 

It is also possible to see from the master equation (9) 
that for an initial Fock state the coherences vanish for all 
times in the Fock-state basis. This is due to the fact that 
the equations of motions for the diagonal elements pk,k it) 
and off-diagonal elements pk,k'(t), with k ^ decouple. 
In particular this means that the equations of motion for 
the off-diagonal elements Pk,k'{t), do not depend on the 
diagonal elements /^^(t). Since for an initial Fock state 
/0fc,fc'(O) = (k\p(0)\k') = 0, the decoupling then implies 
that Pk,k'(t) = for all t and k ^ k' . 

For an initial ground state in the weak coupling 
i(n(t)) <C 1) and for times much shorter than the ther- 
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malization time one obtains 

1 

Po,o = 
Pl,l - 



<"(*)> 
(«(*)>■ 



1 - („(*)) 



(15) 



The analytic solution for the heating function and for 
(n(0)) = is given by 



<n(t)> 



A(t')dt', 



(16) 



where A(i) has been given in Eq. I jlOj l. The exact heating 
function dynamics, given by Eq. fifty, will be used as a 
benchmark for the theoretical quantum simulator results 
in the following subsection. 



with single frequency fields in a pulsed way, the frequen- 
cies covering the relevant part of the environment spec- 
trum of the damped oscillator. We now proceed to show 
how the complete damped harmonic oscillator dynamics 
is recovered when the information about the spectral den- 
sity is used to average the driven oscillator results. We 
first demonstrate that the simulator gives correct heating 
function dynamics and then present the complete solu- 
tion showing that the density matrices of the simulator 
and of the damped harmonic oscillator match. 



1. Heating function 



B. Simulating damped harmonic oscillator with 
driven oscillator 



We have calculated in Sec. |n] the heating function dy- 
namics for a driven oscillator. The oscillator is driven 



Equation © describes the heating function for off- 
resonant driving with ip = 0. When the phase of the 
driving field is taken into account, the heating function 
can be written 



J 



(n){t,u,(p) = \af 



^1 -1 

LUq 



| cos 2 (ip) + cos 2 (u>t + (p) — 2cos(<y9) cos(uit + ip) cos(ujot) — 2— sin(u;t) sin(ajo^) 

I ' ' UJQ 



2 [ — ) sin(w) sin(wi + ip) cos(cj £) + ( — ] |"sin 2 (<p) + sin 2 (ujt + <p)] \. 
\luqJ ' \0J J J 



(17) 



r 



To obtain the final result, we average Eq. (|17fl over the 
phase of the driving field and over the various driving fre- 
quencies by using the spectral density as a weight factor 



7fJ dLJ ^l ^ 7 H< n >(*> w >¥>)-( 18 ) 



(»>(*) 

Averaging over the phase gives, 



<n}(i) 



N 



du)I(u)) 



1 



1+ I — 

UJ Q 



[1 — cos(ujot) cos(wt)] 



uj 



— 2 — sin(wi) sin(wot) 
coo 



(19) 



where N is the normalization factor. We average also 
over the frequency, and thus the normalizing factor, in 
the case we consider, becomes 



2 f°° uj 2 
N = lu c — / duj c 

n Jo oj£ + uj 2 



(20) 



We obtain the result for the heating function dynam- 
ics of a damped oscillator by plugging in the spectral 
distribution [Eq. I|13|) ]. into the phase averaged result, 
Eq. I|19fl . In the regime r < 1 the non-Lindblad-type dy- 
namics dominates. The environment is detuned from the 
system frequency, the most intense part of the environ- 
ment appears for small frequencies uj <C ujq whereas the 
on-resonant part corresponds to the low intensity part of 
the spectral density. We set the parameters as follows: 
r = 0.1, g = 0.045, k B /oj = 80. 

The results from single frequency drives [Figs.n( a )-(c)] 
suggest that the low frequency part of the spectrum (part 
I in Fig. |2J) induces oscillations of the heating function. 
In contrast, the result for oj — oj [Fig. ^ (d)] indicates 
that the near-resonant part of the environment spectrum 
(part II in Fig.|21), produces a Markovian-type heating. 

Indeed, when we calculate the result by integrating 
over the low frequency part, ^ <-j/ujq S* 0.2, the heat- 
ing function displays damped oscillations, see Fig. 0(a). 
When the result is calculated by integrating over the on- 
resonant part, 0.8 ^ t-j/uJa *S 1-2, the heating function 
displays linear dynamics, see Fig. [31 (b) . When we add 
these two results obtained by integrating over the limited 
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FIG. 3: Comparison between the quantum simulator and the exact results for the non-Markovian oscillatory dynamics. The 
results demonstrate how the different parts of the environment spectrum affect the dynamics, (a) < u/u)q < 0.2, (b) 
0.8 < uj/uj < 1.2, (c) 0.0 < uj/luq < 0.2 and 0.8 < uj/lu < 1.2, (d) 0.0 < uj/lu < 1-2. Moreover, the results in (d) show how the 
the random phase of the driving field influences the amplitude of the oscillations of the heating function. Due to the excellent 
match between the simulator and the exact results in (d), it is difficult to see any difference between the two. 



frequency range, the heating function displays the dom- 
inant characteristics of the exact result and the match 
between the two results is already quite close, see Fig. 

(c) . If we then add to the integration the gap between 
the two frequency regimes presented above, the match 
between the quantum simulator and the exact results is 
excellent and it becomes difficult to see any difference 
between the two, see Fig. 0(d). 

These results illustrate how we can identify the origin 
of the dominating features of the heating function dy- 
namics in the various parts of the environment spectrum. 
It is also worth noting that the simulator result in Fig. El 

(d) includes the non-adiabatic effects, which are not con- 
tained in the simple result of Fig. [3] (c). Obviously, the 
inclusion of the non-adiabatic effects enhances the match 
giving excellent agreement between the quantum simu- 
lator and the exact results. Moreover, according to the 
simple displacement scheme for single frequency drives 
(presented in the previous section) , the damping of the 
heating function oscillations arises in a subtle way as an 
indirect consequence of several low frequencies involved. 
All the low frequencies involved affect the simulator dy- 
namics in a similar way when the driving field is switched 
on. However, the effect of the second displacement, when 
the field is switched off, depends on the frequency and it 
consequently damps the oscillations. 

We have also checked the agreement between the exact 
and quantum simulator results for the typical quadratic 
non-Markovian regime (r ~ 10) and for a Markovian 



regime (r > 20, t >> 1/uj c ). For a fundamental theo- 
retical study of various types of dynamics see Ref . |2j| . 

It is also interesting to study the case where the phase 
of the driving field is kept fixed. This illuminates an 
aspect about the effect of the random phase of the en- 
vironmental noise being simulated. The result for the 
fixed-phase, ip = 0, is displayed in Fig. [3] (d), dotted 
line. The random phase does not seem to play a role for 
Markovian heating since the slope of the linear increase 
in (n) for large times is the same for the fixed phase and 
for the phase averaged results. Instead, the phases of 
the low frequency fields affect the short time oscillatory 
behavior of the heating function by influencing the am- 
plitude of the oscillation. Without averaging over the 
phase, and setting ip = 0, gives larger oscillations as can 
be seen in Fig. [3] (d) . When studying the Markovian case 
(r > 20), we noticed, that the fixed phase and random 
phase results very closely agree. This indicates again 
that the random phase plays a significant role only in 
the non-Markovian dynamics. 



2. Density matrix 

So far we have demonstrated that the averaged driven 
oscillator results can mimic correctly the heating function 
dynamics of the open system. It remains to be shown 
that actually the simulated density matrix matches the 
density matrix of the non-Markovian damped harmonic 
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oscillator, given by Eq. (|14fl . 

The single frequency drive (see Eq. Q for the Hamilto- 
nian) creates a coherent state. Following Ref. |22j and us- 
ing the complex field representation of the driving force, 
it can be shown that the amplitude of the coherent state, 
for driving field frequency 10 and phase tp, is given by 



6 ^ 

f3(t,uj,tp) = -nl ^ - - [ 



±L _ 1 

UJ 



\ e iujt — e iuJ ot~l 



where the dimensionless coupling constant k is 

.4 



(21) 



(22) 



The matrix elements in the Fock-state basis are given in 
the usual way 



k-l I QI21 



(23) 



We recall that the analytic result for the damped har- 
monic oscillator is valid for weak couplings. To compare 
the results, we expand the exponential in Eq. I|23|) and 
keep the terms to second order in the dimensionless cou- 
pling constant k. For the sake of simplicity we consider 
here the case of an initial state \n = 0) but the results can 
be generalized in a straightforward way to other initial 
Fock states. Moreover, since we are interested in the non- 
Markovian dynamics, we look for the dynamics for times 
much shorter than the thermalization time. With these 
conditions the diagonal elements of the density matrix, 
up to second order in k, are 



po j0 (t,w,<p) = 1 - \f3{uj,if,t)\ 2 
pi tl (t,w,<p) = \/3(uj,tp,t)\ 2 
pk,k(t,u,<p) = 0, k > 2, 

while the off-diagonal elements are 

p lfi (t,LO,(p) = Pq A = f}(u, <p, t) 
p 2 ,o(t,w,<p) = p^ 2 = (3 2 (uj,ip,t). 



(24) 



(25) 



Similarly to the heating function calculation, the final 
result is obtained by averaging over the phase and fre- 
quency of the driving field. Remembering that I(u>) /N is 
the normalized weight in the averaging and (n)(t, u>, tp) = 
\[3(t,oj,ip)\ 2 , with (n)(t,(j,(p) given by Eq. (|T7jl . we obtain 



1 



1 



Pofl(t) = 1_ T7 / duJ 7T / dcpl(w)\f3(t,w, f)\ 



N 



1 - (n(t)) 



2vr 



2tt 



N 



(n(<)). 



du>- 



2tt Jo 



depI(w)\P(t,u,tp)\' 



Since we know that the heating function of the simulator 
matches the one of the damped harmonic oscillator we 
can conclude that the diagonal elements given by Eq. Ij26|l 
coincide with Eq. (|15fl . In order to prove that the state 
of the simulator is the time dependent thermal state of 
Eq. (|14f> . now we have to prove that the out of diagonal 
elements are equal to zero. 

The phase dependence of pifi(t, w, <p) appears in the 
exponentials exp(±i(p), see Eq. (J2J. Averaging over the 
phase therefore gives 



Pi,o(0 oc 



0. 



(27) 



(26) 



It is worth reminding that Eq. describes the non- 
Markovian dynamics of the damped harmonic oscilla- 
tor in the secular approximation. The same approxi- 
mation in the simulator case amounts at averaging to 
zero the rapidly oscillating terms appearing in p2fl(t), 
see Appendix B. Therefore, in the secular approxima- 
tion, p2.o(t) — 0- Summarizing, we have proved that, in 
the secular approximation, the density matrix of the sim- 
ulator coincides with the density matrix of the damped 
harmonic oscillator. 

It is worth noting that, since the density matrix of the 
damped harmonic oscillator and the one of the simula- 
tor approach coincide, the averaged driven closed system 
mimics correctly also the entropy production typical for 
open quantum systems. In other words, while the closed 
driven oscillator remains always in a pure state (no en- 
tropy production) the averaged state, i.e. the state of 
the simulator, is a mixed state since Trp 2 < 1 for time 
t > and allows the entropy production. To illustrate 
this point further, we display in Fig. 21 by using the simu- 
lator result given by Eq. (|26|l . the von Neumann entropy 
of the system S — —Tr(p\np). The oscillatory behavior 
of the entropy is a direct consequence of the oscillatory 
behavior of the heating function. For a more detailed 
description for the various types of heating function dy- 
namics for engineered reservoirs, see Ref. |23l |. 



IV. IMPLEMENTING THE SIMULATOR WITH 
A DRIVEN TRAPPED ION 

Single trapped ions are a well known example of a real- 
ization of the quantum harmonic oscillator [9( . Recently 
these systems have also been exploited to demonstrate 
the creation of artificial engineered reservoirs for open 
quantum systems 0, ^J. The ability to cool the ion 
to its ground vibrational state and the long coherence 
time make this system ideal for the implementation of 
the quantum simulator described in this paper. 

Important factors for the implementation of the quan- 
tum simulator presented here are the number and fre- 
quency range of the driving fields to be used in the exper- 
iment. Suppose now that we want to simulate the effect 
of non-Markovian off-resonant engineered reservoir to re- 
veal the characteristic quantum mechanical oscillatory 



8 




10 20 30 40 



V 

FIG. 4: The entropy production of the simulated damped os- 
cillator in an engineered reservoir for r = 0.1. The oscillations 
of entropy directly reflect the inherent oscillatory behavior of 
the populations and of the heating function for the short time. 



stands outside the integration and can be taken outside 
of the corresponding summation of the experimental re- 
sults. It is also worth noting that the ambient reservoir 
does not disturb the experiment since its effect occurs on 
few orders of magnitude larger time scale [9| . The vibra- 
tional quantum number of the ion can be measured in 
the standard way . 

The displacement scheme presented in Sec. [H] has 
straightforward implications for interpreting the ion sim- 
ulator results. A sudden switch on and off of the electric 
field displaces the centre of the ion trap due to interac- 
tion between the charge of the ion and the electric field. 
The central point is that in the adiabatic regime, the cen- 
ter of the trap moves slowly compared to the oscillator 
frequency, and the ion follows adiabatically the motion 
of the trap center. This effect plays a key role in the 
damping of the heating function in the experiment. 



V. DISCUSSION AND CONCLUSIONS 



behavior of the oscillator heating function. We set the 
parameters to correspond to Fig. [3 r = 0.1, ui c /2ir = 1.1 
MHz and k B T/io = 80. 

The experimental procedure is the following. We keep 
the step size of the frequency sampling fixed and apply a 
conveniently small number of pulses. We treat the phase 
of the field as a random variable. A typical oscillator 
frequency of the trapped ion in current experiments is 
loq/2-k = 11 MHz [l!j. The frequency range of the pe- 
riodic drive is, e.g., — 12.3 MHz and the separation 
between the frequencies of the sets of drives is 55 kHz. 
In this example the total number of different applied fre- 
quencies is 225. For each frequency, the sets of differ- 
ent duration have to be applied in order to obtain the 
time evolution. The result for the heating function is ob- 
tained by converting the integral of Eq. i|18|) to a sum 
over the sets of applied electric pulses. The result of the 
summation has excellent agreement with the exact and 
theoretical quantum simulator results of Fig. [3] (d) , and 
to avoid repetition we have not plotted in this figure the 
third overlapping curve. Instead it is worth mentioning 
that the x-axis of the figure, for the chosen experimental 
parameters, corresponds to /is time scale. 

The system-reservoir coupling constant is directly pro- 
portional to the intensity of the applied field and the di- 
mensioless coupling constant appearing in the analytical 
solution is directly given by 



where E is the amplitude of the applied field and e the 
charge of the ion. Actually, one can choose the intensity 
for each set of drives ui in the most convenient way for the 
experimental implementation. This is because the used 
intensities can be rescaled when calculating the average 
over the frequencies for the heating function. In other 
words, the intensity dependence \a\ 2 oc E 2 in Eq. i|19|) 



In the past, the average over the stochastic fields has 
been used to describe the effects of ambient reservoirs to 
the oscillator dynamics, e ; g. the heating of the ground 
state of trapped ions 

HE El. Here, we use similar 
displacement operator based formalism to investigate a 
new aspect to the topic. We ask the opposite question: 
how can controlled drives and their average be used to 
model the effects of a reservoir? In fact, we demonstrate 
how to use controlled drives to mimic the effects of both 
the ambient and the artificial engineered reservoirs. 

In addition of the developed quantum simulator, the 
used approach brings to light some interesting features 
of the open system dynamics. We have been dealing 
with a characteristic quan tum mechanical non-Lindblad- 
type master equation |23ll27| . Thus, it is somewhat sur- 
prising that the oscillations of the heating function, a 
typical quantum mechanical dynamical feature, have a 
semiclassical interpretation, when the simulator is im- 
plemented with a trapped ion. Namely, the source 
of the oscillations, according to our interpretation, can 
be seen in the "displacement-coherent state dynamics - 
displacement" scheme described in Sec. [HJ In contrast, 
usually the oscillatory behavior of the heating function 
and non-Lindblad-type-dynamics are connected to the 
appearance of virtual processes [2j| . 

The results also illuminate the role of the random 
phase of the environment fields. In the case of oscilla- 
tory non-Markovian dynamics, the role of the random 
phase is to affect the amplitude of the short time oscil- 
lations. The oscillation is initiated by the low frequency 
fields and its amplitude is given by an average over the 
effects of the random phase of the same frequencies. The 
oscillations are then damped because of the various fre- 
quencies involved. For strictly Markovian case, the phase 
plays hardly any role. 

Our results also shed some light on the problem of the 
existence of a continuos measurement scheme interpre- 
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tation of the quantum trajectories for non-Lindbad-type 
dynamics. In the Lindblad case there is a direct corre- 
spondence between the continuos measurement scheme of 
the environment and the quantum trajectories generated 
by the computer simulations |2l|. For the non-Lindblad 
case a completely satisfactory connection is still lacking 
despite of some attempts |30l 13 1| . and this important 
problem remains open. Thus, it is interesting to note 
that, in order to observe the oscillations of the heating 
function with the quantum simulator described here, the 
act of switching off the field (second displacement) has a 
crucial role. No measurements are allowed while the driv- 
ing field is on. This points to the direction that it is prob- 
able that a satisfactory measurement scheme interpreta- 
tion for non-Lindblad type dynamics does not exist, at 
least not in the framework of conventional quantum the- 
ory. It might be interesting to study in the future what 
the measurement scheme interpretations based on alter- 
native formulations of quantum mechanics 30] would say 
in the context of the quantum simulator proposed here. 

It is important to note that most of the existing ana- 
lytical methods to study open system dynamics rely on 
various approximations, e.g. weak coupling between the 
system and the environment ficj . It would be extremely 
stimulating for the research on open quantum systems to 
possess a general quantum simulator which could be used 
in the regimes where most of the approximations used in 
the analytical calculations fail. We have taken an initial 
step towards this direction and a more general quantum 
simulator for open systems will be a challenging task for 
our future studies. 

In conclusion, we have shown how to use a controlled 
driven harmonic oscillator to simulate the behavior of a 
damped oscillator. In other words, our results demon- 
strate the possibility of studying one of the paradigmatic 
open quantum systems models by means of a closed quan- 
tum system. Moreover, we have discussed in detail the 
implementation of the quantum simulator with single 
trapped ions. In addition, the simulator approach illu- 
minates several interesting aspects of the non-Markovian 
dynamics of a damped oscillator, most notably identify- 
ing the role that the various parts of the non-Markovian 
spectrum play in the system dynamics. 



APPENDIX A 

In this Appendix A we show that for an initial Fock 
state |n = 0) the state of the system as a function of time 
is given by Eq. (|14|) . 

The quantum characteristic function (QCF) of the sys- 
tem at time t can be written |32fl (for a general definition 
of QCF see, e.g., Ref. H) 



(Al) 



where xo is the QCF of the initial state of the system, 
and the quantities Ar(i) and T(t) are defined in terms 
of the diffusion and dissipation coefficients A(t) and ^(t) 
respectively [see master equation ®] as follows 



r(t) 



7(h) dh 



Jo 



(A2) 
(A3) 



Equation (|A1|) shows that the QCF is the product of an 
exponential factor and a transformed initial QCF. 
For a Fock state \n = 0) the initial QCF is |33]| 



xo(0 = e Hm2 . 

By plugging Eq. 1|A4|) into Eq. IjAlf) with £' 
the QCF at time t can be written 

f 



(A4) 



-r(t)/2 



xt(0 = cx p 



A r (t) 



T(t) 



(A5) 



The heating function in turn, for an intial ground state 
of the system, can be written |23l 127] 



Mi)> = ^(e- r(t) -l)+A r (0, (A6) 



which gives 



(n(t)) + i=A r (*) + V r(t) . (A7) 



The right hand side of Eq. (|A7|) appears in the expres- 
sion for the QCF in Eq. (|A5|) which can consequently be 
written 
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Xt(0 = exp 



(n(f)> 



(A8) 



It is also worth mentioning that for high temperature 
reservoir (fcgT/wo 3> 1) the dissipation coefficient ^(t) <C 
A(£) and the heating function dynamics can be approxi- 
mated with Eq. O 0. 

The general expression for a QCF of a thermal state is 

m 



X(£) = exp ^ - 



1 

n+ 2 



(A9) 
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where n is the mean excitation number of the thermal 
state. Comparison between the QCF of the system, 
Eq. QA8JI . and the QCF of the thermal state, Eq. lA9jl . 
shows that the state of the system at each time, for the 
initial state we consider, is a thermal state with chang- 
ing "effective" temperature given by the heating function 
(n(t)). Equation (|14|) follows then directly from this ob- 
servation. 

We note also that the time dependent "effective" tem- 
perature coincides with the temperature of the environ- 
ment only in the long time limit when the system has 
reached a thermal equilibrium with its environment, ft is 
also worth mentioning that the QCF of the thermal state 
with n = 0, Eq. (|A9|I . matches the QCF of the Fock state 
\n = 0), Eq. (|A4(I . since the Fock-state \n = 0) coinsides 
with the T = thermal state. 



We note that 



dip [j3 2 (t,uj,<p)+Pi(t,LJ,ip)] 
dip{c ie - 2 ^ + c 2 e 2lv ) = 0, (B4) 



where ci 2 are functions which do not depend on ip. 
Therefore 



APPENDIX B 

In this Appendix B we show that in the secular ap- 
proximation p2,o(t) — 0. From Eqs. (21) and (25) one 
derives 

P2,o{t) <x- / du I d<pl(ui)[0i(t,uj,<p) + (3 2 (t,uj,(p)} 2 , 
Jo Jo 

(Bl) 



P2,o( f ) x 



du> 



dcpI(u>)Pi (t, lj, <£>)/3 2 (t, w > f) 



DC 
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oc / dujl(uj)- 



■ sin[(tj — uo)t] sin[(w + u!o)t] 







2iuo c \ w c . M , 

1 cos(2o;o<) H sm(2cj £) 

uj J lu q 



(B5) 



with 



Pi(t,UJ,ip) oc 



P 2 (t,u,<p) oc 



LU + OJQ 



lip 



LU — UJq 



( e iujt — e iuJ ot\ 



(B2) 
(B3) 



In the secular approximation, used for the study of the 
damped harmonic oscillator, the rapidly oscillating terms 
appearing in the previous equation average out to zero, 
as one can easily see plotting the real and imaginary part 
of P2,o(0- 
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